On mean field theory for ac-driven elastic interfaces exposed to disorder 
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The analytic description ol ac-driven elastic interfaces in random potentials is desirable because 
the problem is experimentally relevant. This work emphasises on the mean field approximation for 
the problem at zero temperature. We prove that perturbation theory is regular in all orders by 
giving an inductive scheme how to find groups of ill-behaved graphs that mutually cancel, leaving a 
regular expression. In the parameter regimes for which perturbation theory is applicable it agrees 
with numerical results. Further, we determine the dependence of the Fourier coefficients of the mean 
velocity on the parameters of the model. 
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I. INTRODUCTION 

The theory for elastic interfaces in a disordered en- 
vironment, driven by an external dc- force at tempera- 
ture T = is widely understood, and also in the fi- 
nite temperature case some progress has been achieved 
P, i, i, S, i, Hi- At T = the dc driven interface ex- 
hibits an interesting critical point, corresponding to the 
depinning transition. For small constant driving forces 
h, the interface adjusts its configuration to balance the 
driving force and the disorder, but remains pinned and 
does not move on large time scales. If h reaches a critical 
threshold hp, the interface starts to slide with a mean 
velocity v that behaves as v ^ {h~hpf iovh \ hp. The 
critical properties of this non-equilibrium transition have 
been worked out bythe use of functional renormalisation 
group methods 0,11, i, M 

Beyond constant driving forces, experimental achieve- 
ments on the problem of ac-driven elastic domain walls 
in ferroic systems [H, [H, [l^ emphasise the importance 
of a theoretical understanding. Despite the experimental 
progress, and in contrast to the success in understanding 
the dc case, yet there is little advance in the theoretical 
description of the problem, even for T = 0. The exact 
solution of the equation of motion is deemed impossible 
due to the complicated non-linear feedback of the inter- 
face's configuration in the disorder force term. To make 
matters worse, attempts to access the problem for an 
ac-driving force perturbatively by an expansion in the 
disorder strength bring along severe problems [l7j . 

This underlines the importance of the mean-field ap- 
proach, which is the central subject of this article. We 
investigate the behaviour of ac-driven domain walls in 
a disordered environment in the mean-field approxima- 
tion and prove, that the perturbative corrections remain 
bounded in all orders. Further, we indicate, that for large 
enough driving field amplitudes, sufficiently strong elas- 
tic coupling and high frequencies, the perturbative results 
agree very well with the numerics for the full mean field 
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equation of motion. The quantitative statements that 
rely on a special choice for the disorder correlator are 
worked out for elastic manifolds, like for example inter- 
faces between two immiscible fluids or domain walls in 
ferroic systems, exposed to random field disorder. Our 
proof of the regularity of perturbation theory should also 
extend to similar models that describe other interest- 
ing physical systems with disorder [ll, [l^ , for example 
charge-density waves (2^ and fiux lines in type-II super- 
conductors ;2L,i2ai23j]. 

Mean field calculations have been performed for similar 
problems before. D.S. Fisher [13, [111 calculated dynamic 
properties of sliding charge-density waves in a mean-field 
model with dc-driving and argued in favour of a depin- 
ning transition in the strong pinning case. Using smooth 
bounded disorder, he calculated the threshold field for de- 
pinning as well as critical exponents related to the depin- 
ning transition. Furthermore, he considered the response 
in case of an ac-field applied in addition to the dc-driving. 
The perturbation expansion for dc-driven interfaces has 
been investigated by Koplik and Levine |^], who also 
emp hasised on the mean-field problem. Later, Leschorn 
|27| calculated the depinning force and the critical prop- 
erties of the depinning transition for a three-state ran- 
dom field model. Narayan and Fisher Q investigated the 
critical behaviour of charge density waves in the sliding 
regime and worked out the threshold field for scalloped 
disorder potentials. Lyuksyutov considered dynamical 
friction and instability phenomena for the interface mo- 
tion [H. 

The rest of this article is organised in the following way: 
In the next section, we are going to deduce the mean-field 
equation of motion from the original model taken to de- 
scribe the problem of disordered elastic domain walls. In 
section IIIIl we establish the diagrammatic perturbation 
expansion for the mean-field theory and show its regu- 
larity. The bulk part of the inductive proof (the induc- 
tion step) is outsourced to the appendix. After a brief 
review of the problem for a constant driving force in sec- 
tion |lVl we focus our attention to the ac-driving case in 
section |Vl There, we start with sonic qualitative discus- 
sion of the numerical solution and then go on to present 
our attempts to extract information from the first non- 
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vanishing perturbative terms. Due to the comphcated 
non-hnear structure of the expressions involved, numer- 
ical methods had to be employed as well. However, the 
perturbative approach helps a lot to improve numerical 
results. This makes it possible to work out the decay law 
of the Fourier coefficients with their order, in dependence 
of the strength of the driving field. 



II. THE MODEL 

To model I?-dimensional elastic manifolds in a Z) + 1- 
dimensional disordered system, we employ an equation of 
motion that has been introduced in a number of earlier 

works fH,!!!!!! 



the same. Hence, we choose 



dtz{x,t) ^TVlz + h- f{t) + u ■ g{x, z). 



(1) 



The equation does not involve a thermal noise term and 
therefore describes the (classical) system at T = 0. The 
interface profile is described by the single-valued func- 
tion so we do not allow for overhangs. Here, x is 
the i3-dimensional set of coordinates which parameterise 
the interface manifold itself. F denotes the elastic stiff- 
ness constant of the interface, h measures the strength 
of the external driving force and f{t) denotes its time- 
dependence, taken to be of order unity. We did not spec- 
ify f{t) yet to allow for general considerations. The dis- 
order is modeled by g{x,z), its strength is measured by 
u. We assume quenched Gaufiian disorder, characterised 
by its first two cumulants: 



{g{x,z)) = 0, 
{g{x,z)g{x',z')) = {x 



x')A{z^z'). 



(2) 



Here, (. . .) denotes the average over disorder. The func- 
tion A(z) will be specified further down. As the short- 
ranged elastic term suggests, long-range interactions, of 
dipolar type for example, are not covered. 

The corresponding mean field equation (cf. e.g. (25j l 
is obtained via the replacement of the elastic term by 
a uniform long-range coupling. To do this, we have to 
formulate the model ([1]) on a lattice in x-direction, i.e. 
the coordinates that parameterise the interface itself are 
discretised. The lattice Laplacian reads 



D 



^lz{x,) = 



z(xi + acd) + z{xi ~ acd) - 2z{xi) 



D N 

= 51 XI "^y^i [^(^Jd) ~ ^(^')] ' 

where a denotes the lattice constant. To get the mean 
field theory, Jij has to be replaced by a uniform coupling 
but such that the sum over all couplings Jij remains 
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(3) 



Now, the disorder has to be discretised as well, which is 
achieved if we replace the delta function in the correlator 
by (5^(a;, ~ Xj) ^ Jy-a--^/^ (cf. [s^). The result- 
ing equation of motion should be independent of x, just 
the lattice constant a and the dimension enter because 
the disorder scales with a factor a~^/^. Finally, for the 
mean-field equation of motion, we obtain 



dtz^c-[{z)^z] + h- f{t) + i^-g{z), 



(4) 



where c — T /a? and 77 = uja^l"^ . We assume quenched 
Gaufiian disorder with 



(5(z)5(z')> = A(z 



z'). 



(5) 
(6) 



The function A(z — z') shall be smooth, symmetric and 
should decay exponentially on a length scale I. Moreover, 
we require A(0) = 1, as the disorder strength shall be 
measured by 77. For the sake of concreteness, we shall 
choose 



A(z — 2') = exp 



(7) 



whenever we need an explicit expression for calculations. 
This disorder correlator correctly describes the situation 
for an elastic manifold in random field disorder [loj , suf- 
ficiently far away from the critical depinning transition 
point. 

The physical picture of the mean field equation of mo- 
tion is a system of distinct particles, moving in certain 
realisations of the disorder. All of them are harmoni- 
cally coupled to their common mean, i.e. the elastic cou- 
pling between neighbouring wall segments FV^z is now 
replaced by a uniform coupling c • [(z) — z] to the disor- 
der averaged position (z), which in turn is determined 
self-consistently by the single realisations. 

Apart from the correlation length I of the disorder, 
there is another important length scale in the system. In 
the absence of any driving force (i.e. /i = 0), we can easily 
determine the mean deviation of the coordinate z of a 
special realisation from the disorder averaged position 
(z). For = we expect z = 0, at least in the steady 
state and (g]) straightforwardly leads to 

So, r\lc measures the modulus of the average distance 
from the common mean. We will see below, that 77/c is 
an upper bound in the case, that the system moves under 
the influence of a non-zero driving h ^ Q. 

A word on notation: The disorder averaged velocity 
V — (z) will be denoted by the symbol v. 
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III. PERTURBATION THEORY 
A. Diagrammatic expansion 

The differential equation of motion (jj]) is non-linear, 
and due to the influence of the solution on the disorder it 
is impossible to solve it exactly. An ansatz is, to attempt 
an expansion in the disorder strength r]. Therefore, we 
decompose z = Z + where Z = hF{t) (with dtF{t) — 
f(t)) is the solution of the non-disordered problem [rj = 
0) around which we expand, and 



fc=i 



is the perturbative correction. Still, we have the equa- 
tions for Cfc depending on (C)fc, which is also unknown. 
This eventually leads us to a set of two coupled equations 



{dt + c)C = c{0+fj-g{Z + 
9t (0-77- (5(^ + 0), 



(8) 
(9) 



that we can solve iteratively for every order of the per- 
turbation series, if we expand 



9iZ + = 



g(^){Z) 



n=0 



(10) 



If one is interested to keep small orders, this expansion 
of the disorder can only work ii ( i, because £ is the 
typical scale on which g{z) changes. We will come back 



J 



(0, 



to that point, when discussing the special cases for f{t) 
in sections HVl and fVl For the moment, we just do it. 

The propagator corresponding to the left hand side of 
Eq. dH) reads 

G{t) = e{t) ■ e-^K 

Using this propagator, we can formally write down the 
solution and express it order by order in a power series 
in rj. Up to the second order, the solutions are 



(C)i(i) = 0, 



= / dti e 



-c(t-ti) 



(11) 
(12) 



Ci(i) 



t 

(C>2(i)= j dii J dt2e-<'^-'^^A'[Zih)-Zih)], (13) 



t 

Git) = J dtie-^(*-*^) [c (0, ih) + g'{Z{t,)) ■ Ci(ti)] . 



(14) 

Since we assume GauBian disorder, the disorder averaged 
corrections (C)„ vanish for odd n. Due to the nested 
structure, a diagrammatic representation of the pertur- 
bation expansion seems most suited. For the interesting 
quantities (C)fcj the first two non- vanishing orders are 
given by: 



(C)4 






' + 2 





(15) 




• < • • 



r 



The diagrammatic rules are fairly simple: we draw all 
rooted trees with k vertices, and add a stem. Each vertex 
corresponds to a factor g'^™-\Z{t))/m\^ where m counts 
the number of outgoing lines (away from the root). The 
line between two vertices represents a propagator G{t). 
Then Wick's theorem is applied to carry out the disorder 
average. Each two vertices, that are grouped together for 



the average, will be connected by a dashed line. Finally, 
we replace every straight line which, upon removing it, 
makes the whole graph falling apart into two subgraphs, 
by a curly line, corresponding to the propagator of ([9]), 
which is just a Heaviside function 0(t). In the framework 
of equilibrium quantum statistical physics, those graphs 
that involve an internal curly line are called one-particle 
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reducible (IPR). In our classical problem loops only oc- 
cur due to the dashed lines originating from the Gaufiian 
average. 



B. Consistency of the perturbative series 

The perturbation expansion leaves some questions, 
that have to be addressed. It is not immediately obvious, 
that taking the disorder average of p4|) gives the result in 
(fTS]) . i.e. (C)2 {t) = (C2(i))- However, a short calculation, 
using integration by parts reveals this relation to hold. 

Another, much deeper problem is related to the dia- 
grams involving a curly line in their interior. Due to the 
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curly line, they grow linearly in time. In the following, 
we call diagrams non-regular, if they correspond to terms 
which grow unboundedly in time. Koplik and Levine p6l | 
explicitly checked for a time independent driving up to 
sixth order, that the problematic terms of the diagrams 
mutually cancel. We give a very general version, that 
holds for any f{t) and covers all perturbative orders. To 
illustrate, how this works, we present the calculation for 
the fourth order here. The somewhat technical induction 
step, which extends our argument to all orders is given 
in appendix |Al For simplicity, we work with the velocity 
diagrams, that are obtained by just removing the curly 
line from the root. 



2 • 



= / dti e-^(*-*i)A"[Z(i) - Z(ti)] / dt2 / dt3e-=(*^-*^^A'[Z(i2)-Z(t3)] 





tl t2 



dh e-^(*-*^)(-A"[Z(t) - Z{h)]) J dt2 J dts e-^(*^-*^)A'[Z(t2) - ^fe)] 





t tl tl 
S=JdhJ dt2 e-=(*-*^)A"[Z(t) - Z(t2)] J dt3 e-=(*^-*=^)A'[Z(ii)-Z(t3)] 







The modification of the second diagram to express it 
as the sum of the first and S is merely integration by 
parts. The term S now corresponds to the sum of the 
two diagrams. It is not straightforward to see, how S 
behaves generally, but it is easy to see, that it remains 
bounded for large times. Every time integral carries an 
exponential damping term. Basically, we have thereby 
established, that at least up to the fourth order, the per- 
turbation series exists and is well-behaved in the sense, 
that there are no terms that lead to an overall unbounded 
growth in time. 

Our analysis how the cancellations among non-regular 
diagrams (i.e. those that involve an internal curly line) 
generalise to higher orders is presented in appendix 1X1 



that the disorder correlator is cusped, e.g. 



IV. TIME-INDEPENDENT DRIVING FORCE 

Let us consider the mean field equation ^ for the 
special choice f{t) = 1. Particular cases of this problem 
have already been addressed [1, HI US US 113 ■ Assuming 



Ac{z — z') = cxp 



\z ~ z 



(16) 



it is expected, that the system, described by ([U shows 
a depinning transistion. The special feature of a cusped 
disorder potential is, that the resulting disorder force g{z) 
exhibits jumps. At such jumps, the system is pinned 
and a finite threshold force is needed to move it in a 
certain direction. The critical depinning force has been 
determined in Q for the parabolic scalloped potential 
and for a three-state random field model in ^3] ■ In both 
works, the critical exponent was found to be = 1. 

For a disorder correlator given by ([7]), numerical in- 
vestigations suggest, that pinning is exponentially sup- 
pressed for large c. More precisely, the depinning field 
hp obeys hp cx exp[— C • cI/t]] with some numerical pref- 
actor C. This is reasonable, since ry/c is a measure for the 
typical deviation of each single realisation from the mean, 
so for rj/c <ti £, the system is prevented from adopting 
to the minima of the disorder potential by the elastic 
force. Thus, pinning should be diminished. On the other 
hand, for rj/c > £, the equilibrated system for h = 
does place itself at the local disorder minima (for all re- 
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alisations), since on a scale I it is expected to find a local 
minimum. So the existence of a finite perceptable thresh- 
old force is possible for small enough c. Fisher [2^ 
analytically found a threshold field for bounded disorder 
(strong enough) in the case of charge density waves. The 
critical exponent has been found to be /3 = 3/2 for the 
smooth disorder potential, in contrast with /? = 1 for the 
cusped case. Thus, as pointed out in [sj, the exponents 
for the depinning transition in the mean-field case are 
non-universal. 

Well above the depinning threshold, i.e. for h :s> hp, 
perturbation theory should give a good estimate of the 
mean velocity. However, it is not clear, why our pertur- 
bative approach via a systematic expansion in the dis- 
order (essentially equal to the calculations in [2^1) can 
work, since truncating the Taylor expansion (jlOp at a fi- 
nite order is only a good estimate, if ^ <C ^. But, the true 
velocity v is certainly different from h, thence C — {v~h)t 
grows linearly in time. In appendix IB] we illustrate, how 
a resummation of the perturbation scheme leads to a per- 
turbative programme that works. To take only small or- 
ders into account a necessary condition is that 77/c <C ^. 
Choosing ([7]) for the disorder correlator, the first non- 
vanishing order reads 




where we have introduced the function 

4){x) = • X exp (x^) [1 - erf(a;)] (18) 
for convenience. Its asymptotic expansion reads 




For I — > 00, in (|17|) all perturbative corrections vanish, 
hence we get v = h. This was expected, since, if the 
disorder force is correlated over an infinite range, it is 
essentially constant. Taking the average over all possible 
values of the disorder (positive and negative) gives zero, 
hence there is no disorder effect any more. For finite £ the 
velocity is reduced, the smaller c£, the more. For large 
h hp, the difference h — v is expected to be small, 
hence it does not cause much harm if one replaces u by ft. 
on the right hand side of (I17p . rearriving at the explicit 
expansion in the disorder. 

The results from above mainly agree with those in [2^ . 
The major difference is, that Koplik and Levine have 
been working with A^/ (2£) instead of ^ and therefore 
get an expansion in /{£h?), whereas we have a power 
series in rf' /K^ . 
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FIG. 1: Numerical solution of (|19|l for different driving field 
strengths and c = 1.0, r] = 2.5. For the simulation, t and z 
are measured in units such that ui = £ = 1. 



V. CONSIDERATIONS FOR AC DRIVING 
FORCES 

A. Qualitative behaviour and numerical results 

To get an idea about how the system, corresponding 
to the equation of motion with an ac-driving (cf. (|^) 

dtz = c ■ [{z) — z] + h ■ cosLvt + 7] ■ g{z) (19) 

behaves, we implemented a numerical approach. The 
disorder is modelled by concatenated straight lines, the 
values of the junction points are chosen randomly from 
a bounded interval. The correlator has been checked to 
be perfectly in agreement with ([7]). 

Before discussing the numerical trajectories, we note a 
first property of the equation of motion (jl9p . It contains 
a symmetry of the (disorder averaged) system, namely 
that all disorder averaged quantities are invariant under 
the transformation h —h and z ^ —z, which implies 
V —V. We have hereby fixed the initial condition to 
be z(0) = for all realisations. If one chooses another 
initial condition, its sign has to be inverted as well, of 
course. In the steady state, i.e. for t :§> c^^ (as c^^ is the 
time-scale on which transience effects are diminished, see 
below), the trajectory must therefore obey the symmetry 
h — !■ —h, V — !■ —V. This symmetry is obviously reflected 
in the numerical solutions (see fig. [1]). 

An interesting consequence of this symmetry is, that 
the even Fourier coefficients of the solution v{t) (which is 
periodic with period 27r/w) vanish. Once the steady state 
is reached, the symmetry requires v{t) = —v(t + tt/uj). 
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2: Numerical solution of equation 
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for h = 6.0 



c = 1.0 and r\ — 2.5 for different frequencies, t and z being 
measured in units such that ljq = £ = 1. 




FIG. 3: Numerical solution of equation (|19|l for different elas- 
tic constants and h = 40.0, 77 = 10.0. The units of t and 2 are 
chosen such that ui = £ = 1. 



For the even Fourier modes this means 



C2N = j dit;(i)e'2^"* 


JL 2L 

= Jdt z;(t)e'2^"* + Jdtv{t+^) e'2^"* = 0. 


The typical picture of a v-h-plot is that of a single hys- 
teresis for /i ^ 77 and a double hysteresis for /i ^ 77. In 
an intermediate range, we find a single hysteresis with a 
cusped endpoint. The qualitative shape of the solution 
trajectories agrees with numerical results [sil [33 |. that 
have been obtained as solutions for ([1]) in the case of finite 
interfaces with periodic boundary conditions. Moreover, 
as the frequency is sent to zero w — > 0, the hysteretic tra- 
jectory approaches the depinning curve for an adiabatic 
change of the driving field. This is shown in fig. O 

In the following, we want to give a qualitative discus- 
sion of the hystereses in the case of small elasticity c. 

Weak fields h <^ rj: In the case of weak driving fields, 
the typical system in a certain disorder configuration re- 
mains in a potential well of the disorder. The elastic force 
may slightly shift the centre around which z{t) oscillates, 
but this is not very important, since we can instead think 
of an effective potential. To understand the hysteretic be- 
haviour, it is instructive to think of the force field g(z) 
instead of the potential. Starting at h{t) — for large 
enough t (i.e. in the steady state) we expect a certain 
realisation to be located at the zero point 5(2:0) = of a 
falling edge, since this corresponds to a stable configura- 
tion. As the field grows, the system starts to move in the 
direction of growing z, where the disorder force competes 
with the driving. Because in the vicinity of the poten- 
tial minimum, the disorder force g{z) behaves approxi- 
mately linear in z, the acceleration is approximately zero 
and the velocity almost constant. This changes when the 
driving is about to reaching its maximum. The slower 



the growth, the smaller the velocity. At the maximum, 
the velocity equals zero, as the driving and the restitu- 
tional disorder force compensate. For decreasing h{t), 
the restitution force wins and pushes the system back 
in the direction of the potential minimum. Hence, the 
velocity v turns negative short-time after the field has 
reached its maximum and is still positive. Once the sta- 
ble position zq is reached again, the same starts in the 
negative direction. 

Certainly, the restitutional disorder force need not con- 
tinuously grow with z, but may exhibit bumps or similar 
noisy structure, but those details average out when tak- 
ing the mean over all disorder configurations. 

Strong fields h ^ rj: In the case of strong driving 
amplitudes, we encounter the situation of a double hys- 
teresis. Again, starting at h(t) = for i ^ c~^, we as- 
sume the system to be located at the zero g{zQ) = of a 
falling edge of the (effective) disorder force field. As h{t) 
grows, we first have the same situation as in the case of 
weak driving: The disorder acts restitutionally and thus 
keeps the velocity small and leads to a small slope dv/dh. 
Once the field is of the order 77, the typical maximum of 
a disorder force, the system is no longer locked into a 
potential well, but a cross-over to sliding behaviour sets 
in. On further increasing h, the system finally arrives 
at a slope s = dv/dh, which depends mainly on 77 and 
c. After the field reaches its maximum, the velocity de- 
creases with the field, the slope being dv/dh « s, also 
if this slope has been different just before the field am- 
plitude has been reached. This slope approximately re- 
mains, until the field is weaker than the typical disorder 
force, when the system is again trapped in a potential 
well. Since on rising edges of the disorder force, driv- 
ing and disorder point in the same direction, the system 
will rarely sit there (it moves away very fast). The ve- 
locity becomes negative before h — 0, since the system 
shdes down the faUing edge {dg/dz < 0) of the disorder 
force. At /i = everything starts again in the negative 
direction. An example for fairly large field amplitudes is 
shown in fig. [31 
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So far, our discussion has emphasised on small c by 
absorbing its effect into an effective disorder picture. The 
effect of larger c is to couple the configuration z{t) of 
every realisation strongly to the mean {z{t)). This wipes 
out the effect of disorder in the time regime when h{t) 
takes on small values. As we have discussed above, in 
those time intervals the possibility to explore the shape 
of the individual disorder landscape plays an important 
role. Thus for larger c the double hysteresis winds around 
a straight line, connecting the extremal velocities. This 
can be seen in fig. [H 
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B. Validity of perturbation theory 

For an oscillating driving force, the question is still 
open, whether one may assume C to be small compared 
to £. If c is large, any particle moving in a particular 
realisation of a disorder potential is strongly bound to 
the disorder averaged position. This prevents it from 
exploring the own disorder environment and thus large 
c effectively scale down 77. All realisations stay close to 
the disorder averaged position, the mean deviation being 
rj/c. A problem now occurs, if the disorder averaged 
position deviates strongly from the 77 = solution. For 
/i 3> 77 this can only happen during those periods, where 
h cos ujt takes on small values. The time, that has to 
elapse, until every system has adopted to its own disorder 
realisation, and hence the time until the system can be 
pinned, is (see below). For perturbation theory to 
work, this time must be large compared to the length 
of the period during which h < r], which we roughly 
estimate as r]/(ioh). This gives us a second condition for 
the applicability of perturbation theory: h/rj ^ c/uj. 

In summary, the conditions for perturbation theory to 
hold are the following. The driving force amplitude h 
has to be large compared to 77, h/rj ^ maxjc/u;, 1} to 
make the series expansion work and to guarantee that 
the disorder averaged solution stays close to the 77 = 
trajectory (around which we expand). Moreover, c must 
be large (c 3> Tj/f,) to ensure proximity of each realisation 
to the disorder average. 



C. Simple perturbative estimates 

Before embarking on conclusions that can be drawn 
from the perturbative solution of the equation of motion 
(HI, we shall determine the typical deviation of the po- 
sition of a single realisation z{t) from the mean {z{t)). 



FIG. 4: Comparison of the numerical solution of equation (|19|) 
with the result obtained from the first non-vanishing pertur- 
bative order for (1) h = 3.0, c = 3.0, 77 — 1.5 and (2) h = 1.0, 
c = 1.0, ?7 = 0.6. The units of t and z are chosen such that 

UJ=t=l. 

This has the following bound 

{{{z) - zf) ^ r^^Cl) + 0{r,') 
t 

^rf j diidt2 e-^(2t-ti-t2)^[^(^^) _ z{t2)] 


t 

< J dhdh e-^(2t-ti-t.) = J (1 - e-=*)' . 


(20) 

The estimate simply replaces the disorder correlator A 
by its maximum and therefore gives an upper bound. 
Strictly, it is only true to order 0{rf'). This confirmes 
our claim from section |TT1 that an estimate for the mean 
deviation from the averaged solution is given by rj/c. 

D. Perturbative harmonic expansion 

As has been discussed in section IV B[ for h/r] ^ 
maxjc/cj, 1} and t^/{c1) <^ 1, perturbation theory should 
do pretty well. A direct comparison, shown in fig. 31 
confirms an excellent agreement. So, at least for weak 
disorder, when perturbation theory is valid, one can hope 
to extract some information from the lowest order. For 
an ac driving force, even this lowest perturbative order 
for the velocity is a very complicated expression. The 
diagrammatic prescription yields up to the order 0{rf') 

t 

v{t)^hcos.ujt + Tf j di'e-"(*-*')A'[Z(t)-Z(i')]. (21) 


Remember, that Z{t) — (h/uj) smut is the solution for 
the problem without disorder, around which we expand. 
It seems reasonable to aim a harmonic expansion of the 
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mean velocity v. The ansatz therefore is 

oo 

v{t) = ^ [gn COS Nujt + bN sin Nivt], (22) 

Ar=l 

for N odd. Recall, that, for reasons of the h —h 
and V — !■ —V symmetry of the trajectory, which has been 
discussed in the previous section, the Fourier coefficients 
for even N vanish. 

Starting from the first order result for v (|21|) . we ex- 
press the disorder correlator by its Fourier transform 

A'[Z{t) ^ Z{t')]= f ^(ig)A(g)e'«S[^'"'^*-"'"'^*'l (23) 
J 27r 

and expand the exponential term in a double Fourier se- 



ries in t and t', respectively: 



Aa sin ut 



n— — oo 



^-cit-t')-iasinu>t' ^ Jni-a)- 



c + into 



Here, Jn(a) are the Bessel functions of the first kind. As 
we are interested only in the behaviour for large enough 
times (the steady state solution), we remove all terms 
that are damped out exponentially for t ^ from the 
very beginning. Note, that is indeed the time scale 
for the transience, as has been claimed before. 
For the mean velocity, we obtain 



v{t) ^ hcosojt + T]'^ 22 / 2^(^^)^(^)'^™ ( ) ( ) ^ 

m,n— — CX3 ^ ^ ^ ^ 

I 



In principle, this is already a Fourier series representa- 
tion, not very elegant, though. The argument (111 + 12)^1 
of the expansion basis exponentials promises a rather 
complicated structure for the coefficients. A first observa- 
tion, however, can already be made: Under the q integral 
we find an odd function (ig)A(q) and a product of two 
Bessel functions of order m and n, respectively. For the 
q-integral to result in a finite value, a function is required 



that is not odd in q. This necessitates the product of the 
two Bessel functions to be odd, or, equivalently, m + ri 
to be an odd number. Whence, we conclude, that to first 
perturbative order, our symmetry argument (Fourier co- 
efficients for even N must vanish) is fulfilled exactly. 

It requires some tedious algebra to collect all contri- 
butions belonging to a certain harmonic order from the 
double series. Eventually, we obtain the series expansion 



v{t) 



— — cos ujt + ^ [Ajsi cos Nujt + Bm sin Nujt] 
^ v=i 

E(-l)"ni4'jv-„,„ nKN+ n,n 

1 + n^uP-lc^ ^ 1 + n^bj^lc^ 



A 



N 



if Up' 



77^ LO 



ri=l 



Kn-N,n 

n=\ ' n=0 ' n=N ' 



St 



K 



N+n,n 



00 

St 



^ n—i\ ,71 



(25) 
(26) 
(27) 



For convenience, we introduced the following abbreviation (depending only on the parameter ratio h/{uj£)), in which 
3-^3 denotes the generalised hypergeometric function [33]. 



=2 /^M'"^"^' r([m + n + 2]/2) ^^ 



3F3 



^/t: • to! • n! 
TO + n+1 TO-|-n-i-2 TO + n + 2 



(28) 
(29) 



, {to + l,n+l,TO + n+ l},- 



2p 



Note, that taking — > is forbidden here, as we used uj ^ while deriving the coefficients and moreover per- 



9 



/i = 2.7; c = 1.0; r; = 0.6 




10 20 30 40 50 

N 



FIG. 5: Plotting the logarithms of \An\ and \Bn\ reveals the 
exponential decay with A^. In the regime where numerical 
errors do not dominate the result, a linear regression seems 
appropriate. 

c= 1.0 ; J/ = 0.6 



a H 




J I \ \ I I \ L 

1 2 345678 

h 



FIG. 6: Performing the linear regression for many h yields 
slopes Q and /3 appearing to depend on /i in a power-law 
fashion. 

turbation theory breaks down (recall that h/r] 3> c/u). 
The same holds for £ —^ 0. The remaining extreme limits 
cij — !■ cxD and £ ~^ oo are not interesting, since in these 
limits the disorder is rendered unimportant. Therefore, 
in the following, we assume finite (positive) values for £ 
and Lv and moreover set them equal to one lu — £ — 1, by 
appropriately choosing the units for z and t. 

Now, we are left with three dimensionless parameters: 
h, c and 77. The dependence of the first order perturbative 
Fourier coefficients on 77 is trivial. The dependence on c 
is also evident, as can be read off from (|26l27p . For larger 
c, the system is more tightly bound to the non-disordered 
solution, supressing perturbative corrections. 

The most interesting but also the most difficult is the 
dependence of the Fourier coefficients on h. Actually, 
there are two competing effects. On the one hand, large 
driving strengths render the disorder unimportant in all 
cases accessible through perturbative methods. In a nut- 
shell, the expansion parameter is r]/h. On the other 
hand, if one thinks of g{Z(t)) as a function of time, the 
more rapid Z{t) changes the more g fluctuates on short 



time scales and thus brings higher frequency contribu- 
tions to v{t). The first remark is reflected in the overall 
weight of the Fourier coefficients as corrections to the 
non-disordered case, decreasing with h. The second idea 
is expected to express itself in the decay of the Fourier 
coefficients with N. The larger h, the weaker we expect 
this decay to be. 

In equations pS)) and (pS)) . the dependence of the 
higher harmonics on h is hidden in the Km,n as func- 
tions of h: The maximum of the Km^n as functions of 
the parameter h/(uj£) shifts to larger values as m or n 
increase. In the integral representation this can be seen, 
as the Bessel functions take their first extrema at large 
arguments for large indices. However, the complicated 
way in which the Kra,n functions enter Ajq and _Bjv hin- 
ders an analytic access to the decay law. A numerical 
determination of the Fourier coefficients for the pertur- 
bative result reveals an exponential decay, cf. fig. [5] 
The noisy behaviour for > 40 is due to numerical fluc- 
tuations. Note, that these fluctuations are of the order 
10~^^, which is quite reasonable. The plot in fig. [5] is 
mere illustration of a more general phenomenon. This 
exponential decay has been found for many sets of pa- 
rameters, thus one is led to the ansatz 

|^jv|~e~"^ ; |i?w|~e-'5^, (30) 

where a and (3 can be estimated through a linear regres- 
sion up to a suitable A'max- Of course, it is not expected, 
that a and (3 are distinct, nor that they depend on the 
parameters in different ways. Determining both just dou- 
bles the amount of available data. 

As our results are first-order perturbative, a and (3 
must not depend on 77. The main interest now focusses 
on the dependence of the decay constants on h. The 
results from a linear regression for a series of /i- values, c 
and rj kept fixed, suggest a power-law dependence 

a{h) ^ ■ h'^" , (3{h) = Cp ■ h'^^ . (31) 

Fig. [S] displays this relation for a particular example. 
Repeating this data collection and subsequent regression 
for different values for c and rj yields the results sum- 
marised in table [J While the exponent ^ appears con- 
stant ^ K, 0.6, the prefactor seems to depend on c. An 
attempt to redo the same procedure, done for h, with 
the parameter c to gain information about the functional 
dependence of a and /3 on c yields a complicated but 
rather weak dependence, which gives no further insight. 
The obvious approach is to visualise the dependence of 
Ca and Cfj on c. The linear fit in fig. [7] gives a fairly 
tiny slope, so the dependence of the decay constants on 
c may be assumed to be weak. 

Certainly, it is desirable to ascertain the validity of this 
decay law beyond perturbation theory. In a few words, 
it ought to be explained, why we have not been able to 
do it. First of all, the logarithmic plots of the Fourier co- 
efficients in fig. O exhibit fluctuations around the linear 
decrease. This "noise" is authentic and not attributed 
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TABLE I: Results for the regression ((3T|) . 



c 


V 










1.0 


0.6 


1.52 


1.52 


0.61 


0.61 


1.5 


0.6 


1.56 


1.56 


0.58 


0.59 


2.0 


0.6 


1.56 


1.59 


0.58 


0.60 


2.5 


0.6 


1.63 


1.62 


0.61 


0.61 


3.0 


1.0 


1.66 


1.66 


0.63 


0.63 


o.o 


l.U 


1 71 
1.(1 


l.Do 


U.DZ 


U.DZ 


4.0 


1.0 


1.69 


1.65 


0.61 


0.60 


4.5 


1.0 


1.72 


1.68 


0.62 


0.62 


5.0 


2.0 


1.71 


1.67 


0.61 


0.60 


5.5 


2.0 


1.73 


1.73 


0.61 


0.62 


6.0 


2.0 


1.77 


1.72 


0.62 


0.62 


6.5 


3.0 


1.76 


1.77 


0.62 


0.62 




c 



FIG. 7: Plot of the change of the prefactor C(c) in pip on c. 
The linear fit yields a fairly tiny slope. 

to numerical inaccuracies. The exponential decay of the 
Fourier coefficients is superimposed on a true, compli- 
cated dependence. Hence, it requires a lot of data points 
to obtain reasonable data. Since the Fourier coefficients 
for even N vanish, in the example of fig. [5] the regression 
can be carried out over around 15-20 data points. This is 
a fair number. The quality relies heavily on the accuracy 
of the numerical determination of the Fourier coefficients. 
In Fourier analyses of the numerics for the full equation 
of motion (jl9p , we did not manage to get a precision bet- 
ter than of the order of 10^'^. This means, the regression 
has to be stopped at A^max, where logAN,„^^ « —7. In 
the example of fig. [5l this leaves us with less than 5 
data points. In view of the natural fluctuations, a linear 
regression is not sensible any more. 



VI. CONCLUSIONS 

The mean- field version of the problem of ac-driven elas- 
tic interfaces in disordered media admits a regular per- 
turbative treatment that, where applicable, agrees very 
well with the numerics for the full equation of motion. It 
has been shown, that diagrammatic contributions with 



an unbounded increase in time cancel among each other, 
leaving a well-behaved perturbative expansion. 

The solutions to the mean-field equation of motion 
are found to share many features with the numerical so- 
lutions of the original problem, like the hysteretic be- 
haviour of the u-/vplots. 

Unfortunately, the perturbative expressions are very 
complicated and thus only of little use for analytic in- 
sights. However, they improve numerical results tremen- 
dously which made possible to establish the dependence 
of the decay constants of the Fourier modes on /i as a 
power law a,(3{h) = C{c) ■ h^^ with ^ ~ 0.6. 
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APPENDIX A: REGULARITY OF THE 
PERTURBATION EXPANSION 

In section UlI Bl we have analysed, how the unbounded 
contributions, contained in the two diagrams that involve 
a curly line, mutually cancel in the second non- vanishing 
perturbative order. In this appendix, we are going to 
explain how this cancellation process generalises to all 
orders in perturbation theory. As before, for simplicity, 
we work with the diagrams for the disorder-averaged ve- 
locity, that arise by just removing the curly lines from 
the root of the diagrams for (C) (cf. equation ^T5\i ). In 
a velocity diagram contributing to the n-th order (recall, 
that only for even n the corrections are non-zero), any 
curly line connects two trees of order p and q (both even) 
with the restriction p + q = n. Both trees appear in the 
expansion of lower orders, namely p and q, respectively. 
In the following, we want to sketch an inductive proof 
for the claim that the unbounded terms originating from 
trees with curly internal lines cancel among each other. 

Let us assume, that for order n we have achieved to 
ensure regularity. For every unbounded tree T, there 
is thus a set T^,...,T° of, let us call them cancelling 
trees, such that T + + . . . + T°- is a. regular, bounded 
expression in time. As a starting point for the induction, 
take n = 4, where the validity of the claim has been 
verified in section llll Bl It is now the task to validate the 
regularity for order n + 2. First of all, we consider the 
process of attaching the root of a regular tree S (with no 
internal curly line) of order s by a curly line to a vertex v 
of another regular tree R of order r = n + 2 — sto obtain 
a new irregular tree A of order n + 2. The vertex v must 
be connected to another vertex w € R hy a dashed line, 
to carry out the GauBian disorder average. Without loss 
of generality, we assume that v is connected to w by a 
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path that first makes a step towards the root. The rules 
for the diagrammatic expansion ensure, that there is a 
maximal regular subtree T d R, which contains v and w. 

Using partial integration, it is possible to move the 
vertex to which S is connected (via the curly line) to a 
neighbouring vertex in T. Thus, it is possible to move 
the connection vertex along the unique way (in T) from 
V to w. We are going to show, that once w is reached, 
we have obtained the cancelling tree which is unique. 
Diagrammatcially, the process of moving the connection 
vertex from v to w reads: 



D 



Here, the blank circle represents S, the lightgrey circle 



J 



stands for the subtree i?i of R, to which v connects and 
the darkgrey shaded circle denotes trees which run out 
of V (summarised in the following as i?2)- Certainly, in 
general there may be dashed lines between the dark- and 
the lightgrey circle, which we have omitted as they are 
not relevant for the forthcoming discussion. The dotted 
line just serves as a joker - it is not important to specify 
how many trees go out of v. The last term D collects the 
left-over terms from the partial integration. Note, that, 
if it takes several steps to go from v to w, the interme- 
diate expressions (in the partial integration) are no valid 
diagrams. 



To illustrate the procedure, we take a look at the first 
step: 



o- 



Ti 



R,{t) dhe~<^^-'^\-irA('^+''^[Z{T)-Z{h)]R2{h) UhS{h) 



Ti 



Ri{t) / dt2S{t2) / dtie-=(^^-*^)(-l)'^A(^+'')[Z(T) - Z{h)]R2{h) 



Ti ti 

-Ri{t) J dh e-^(^i-*i)5(ti) J dt2e-^(*i-*^)(-l)''A(^+^)[Z(T)-Z(i2)]i?2(i2) 



I 



The order of the derivative (i.e. the number of outgoing 
lines) of w and v are denoted by /i and i/, respectively. 
The time, at which the whole diagram is to be evaluated, 
is the time corresponding to the vertex to which v is 
connected is given by Ti, ti is thus the time associated 
to V and so on. The time of w is t. Thus, we see, that 
if w is not the vertex to which v is directly connected 
(then Ti ^ T in general), the first expression after par- 
tial integration cannot be a valid diagram: v has lost one 
order of derivative (z/ — 1 lines go out instead of v) , but 
the derivative of the correlator A has not changed. A 
valid diagram is then obtained, when the connection of 
S has reached w. Then, v has lost an outgoing line, but 
w received one more and we indeed have achieved a can- 
celling tree: the factor (—1)'' remains, the true diagram, 
however, has (— 1)"^"-^. The signs are different, thus the 
two trees cancel. The left-over term from the partial in- 
tegration is again regular, as can be seen since all time 
integrals carry an exponential damping term. It is clear, 
that this is generally true for every partial integration 
step. 

To go one step further, we assume now S to be irregu- 
lar. Essentially, the same procedure works, but there are 
more cancelling trees: one has take all cancelling trees 
{S^} for S into account (which exist by induction hy- 



pothesis), thus S is replaced by J2^^ ^^'^ thence the 
left-over terms are again regular. 

A possible irregularity of R can be accounted for in 
the same way. It is, however, important to explain why 
this is possible, i.e. what is v and w in the cancelling 
trees for R. In the case of irregular S the problem was 
easy, since all trees have a unique root. As we have seen 
already, the procedure of creating cancelling trees does 
not change the structure of regular subtrees. Thence, all 
cancelling trees for R contain T. This makes clear, which 
V and w have to be chosen in the cancelling trees: they are 
well-defined in T and T is a well-defined subtree of the 
cancelling trees. Thus, repeating the whole procedure 
described above for all cancelling trees of R yields the 
complete set of cancelling trees for A in the most general 
setting. 



APPENDIX B: RESUMMATION OF THE 
PERTURBATION EXPANSION 

For a constant driving force, it is not clear, why a 
perturbative approach using (jlOp should work. So, we 
approach the perturbation expansion from another di- 
rection, which will turn out to be a resummation of our 
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former expansion in the disorder. Take our original equa- 
tion of motion 

dtz^c- [{z) -z] + h- fit) + 77 • g{z) (Bl) 

and decompose z = X + where X = (z). Thus (^) = 0. 
This gives us two non-hnearly coupled non-linear differ- 
ential equations 

dtX{t)^h-f{t)+rj-{g{X + 0) (B2) 
{dt + c)m =v[9iX + 0-{9iX + 0)] (B3) 

If c is large enough, one can always achieve ^ ^ ^ and 
thus a Taylor expansion of the disorder around X{t) keep- 



ing only lowest order-terms seems reasonable. Instead of 
a systematic expansion in the disorder, we now perform 
power-counting in ^. This leads to a recursive structure 
for 



g{X) + g'{X)i - {g'{X)i) + -g" {X)e + 



and thus a self-consistent equation for X{t) 



J 



00 ^ 

aiX = /i./(0 + 7y^-(g(")(X)r 



= h ■ fit) + ii{g{X) + w'(^) / dt'e~^(*-*') [g{X) + ...]+...) 



The graphical structure is now similar to that of section UlI Al It reads 



dtX =h- f{t)+>i^^ -H3 





This series only consists of "irreducible" (IPI) graphs, 
where no line can be cut such that the whole graph falls 
apart into two and no curly lines occur. Otherwise the di- 
agrammatic rules are essentially the same as before. Ev- 
ery vertex corresponds to (7^™'*(X)/m!, where m counts 
the number of outgoing lines. The full lines are propaga- 
tors of the differential equation (jB3|) for ^. 

On analysing the expansion of the first and simplest 
diagram in the disorder (up to Oirf)) 



X ^ = r +2 



g{X) = g{Z + (C» = g{Z) + g'{Z) (C> + . • . , 

one inspects that the "new" diagrammatic expansion 
(|B4I) is merely a resummation of our old one. However, if 
a constant driving force is exerted, ^(i) remains bounded 
(and, depending on c, also small) at all times, in con- 
trast to C(^)- This can be seen by confirming the earlier 



estimate for the mean deviation of a realisation from the 
mean position (cf. section [TT|, which is exactly S^. We 
observe that 



{{{z) ^ z)') = {e) 



Oiv% (B5) 



As < 1, to this order 77/c is an upper bound for the 
typical distance from the mean. Thus this perturbation 
expansion should work fine for ry/c <C ^. As far as there 
is overlap, this also agrees with the result in [2^. For a 
constant driving force h, the first order correction to the 
velocity yields a self-consistent integral equation 



v^h + f]^ dt' e-^(*-*') A'[w(t - t')] + 0{r]^), (B6) 



which has been used in computing equation (|17p . 
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